Antibody production relies on the tRNA inosine wobble modification to meet biased codon demand.

Antibodies are produced at high rates to provide immunoprotection, which puts pressure on the B cell translational machinery. Here, we identified a pattern of codon usage conserved across antibody genes. One feature thereof is the hyperutilization of codons that lack genome-encoded Watson-Crick transfer RNAs (tRNAs), instead relying on the posttranscriptional tRNA modification inosine (I34), which expands the decoding capacity of specific tRNAs through wobbling. Antibody-secreting cells had increased I34 levels and were more reliant on I34 for protein production than naïve B cells. Furthermore, antibody I34-dependent codon usage may influence B cell passage through regulatory checkpoints. Our work elucidates the interface between the tRNA pool and protein production in the immune system and has implications for the design and selection of antibodies for vaccines and therapeutics.

Inherent genome instability underlies trisomy 21-associated myeloid malignancies.

Constitutional trisomy 21 (T21) is a state of aneuploidy associated with high incidence of childhood acute myeloid leukemia (AML). T21-associated AML is preceded by transient abnormal myelopoiesis (TAM), which is triggered by truncating mutations in GATA1 generating a short GATA1 isoform (GATA1s). T21-associated AML emerges due to secondary mutations in hematopoietic clones bearing GATA1s. Since aneuploidy generally impairs cellular fitness, the paradoxically elevated risk of myeloid malignancy in T21 is not fully understood. We hypothesized that individuals with T21 bear inherent genome instability in hematopoietic lineages that promotes leukemogenic mutations driving the genesis of TAM and AML. We found that individuals with T21 show increased chromosomal copy number variations (CNVs) compared to euploid individuals, suggesting that genome instability could be underlying predisposition to TAM and AML. Acquisition of GATA1s enforces myeloid skewing and maintenance of the hematopoietic progenitor state independently of T21; however, GATA1s in T21 hematopoietic progenitor cells (HPCs) further augments genome instability. Increased dosage of the chromosome 21 (chr21) gene DYRK1A impairs homology-directed DNA repair as a mechanism of elevated mutagenesis. These results posit a model wherein inherent genome instability in T21 drives myeloid malignancy in concert with GATA1s mutations.

Myeloid cell subsets that express latency-associated peptide promote cancer growth by modulating T cells.

Myeloid suppressor cells promote tumor growth by a variety of mechanisms which are not fully characterized. We identified myeloid cells (MCs) expressing the latency-associated peptide (LAP) of TGF-ß on their surface and LAPHi MCs that stimulate Foxp3+ Tregs while inhibiting effector T cell proliferation and function. Blocking TGF-ß inhibits the tolerogenic ability of LAPHi MCs. Furthermore, adoptive transfer of LAPHi MCs promotes Treg accumulation and tumor growth in vivo. Conversely, anti-LAP antibody, which reduces LAPHi MCs, slows cancer progression. Single-cell RNA-Seq analysis on tumor-derived immune cells revealed LAPHi dominated cell subsets with distinct immunosuppressive signatures, including those with high levels of MHCII and PD-L1 genes. Analogous to mice, LAP is expressed on myeloid suppressor cells in humans, and these cells are increased in glioma patients. Thus, our results identify a previously unknown function by which LAPHi MCs promote tumor growth and offer therapeutic intervention to target these cells in cancer.

Clonal selection of stable aneuploidies in progenitor cells drives high-prevalence tumorigenesis.

Chromosome gains and losses are a frequent feature of human cancers. However, how these aberrations can outweigh the detrimental effects of aneuploidy remains unclear. An initial comparison of existing chromosomal instability (CIN) mouse models suggests that aneuploidy accumulates to low levels in these animals. We therefore developed a novel mouse model that enables unprecedented levels of chromosome missegregation in the adult animal. At the earliest stages of T-cell development, cells with random chromosome gains and/or losses are selected against, but CIN eventually results in the expansion of progenitors with clonal chromosomal imbalances. Clonal selection leads to the development of T-cell lymphomas with stereotypic karyotypes in which chromosome 15, containing the Myc oncogene, is gained with high prevalence. Expressing human MYC from chromosome 6 (MYCChr6) is sufficient to change the karyotype of these lymphomas to include universal chromosome 6 gains. Interestingly, while chromosome 15 is still gained in MYCChr6 tumors after genetic ablation of the endogenous Myc locus, this chromosome is not efficiently gained after deletion of one copy of Rad21, suggesting a synergistic effect of both MYC and RAD21 in driving chromosome 15 gains. Our results show that the initial detrimental effects of random missegregation are outbalanced by clonal selection, which is dictated by the chromosomal location and nature of certain genes and is sufficient to drive cancer with high prevalence.

Quantitative mapping of the cellular small RNA landscape with AQRNA-seq.

Current next-generation RNA-sequencing (RNA-seq) methods do not provide accurate quantification of small RNAs within a sample, due to sequence-dependent biases in capture, ligation and amplification during library preparation. We present a method, absolute quantification RNA-sequencing (AQRNA-seq), that minimizes biases and provides a direct, linear correlation between sequencing read count and copy number for all small RNAs in a sample. Library preparation and data processing were optimized and validated using a 963-member microRNA reference library, oligonucleotide standards of varying length, and RNA blots. Application of AQRNA-seq to a panel of human cancer cells revealed >800 detectable miRNAs that varied during cancer progression, while application to bacterial transfer RNA pools, with the challenges of secondary structure and abundant modifications, revealed 80-fold variation in tRNA isoacceptor levels, stress-induced site-specific tRNA fragmentation, quantitative modification maps, and evidence for stress-induced, tRNA-driven, codon-biased translation. AQRNA-seq thus provides a versatile means to quantitatively map the small RNA landscape in cells.

RAD21 is a driver of chromosome 8 gain in Ewing sarcoma to mitigate replication stress.

Aneuploidy, defined as whole-chromosome gain or loss, causes cellular stress but, paradoxically, is a frequent occurrence in cancers. Here, we investigate why ∼50% of Ewing sarcomas, driven by the EWS-FLI1 fusion oncogene, harbor chromosome 8 gains. Expression of the EWS-FLI1 fusion in primary cells causes replication stress that can result in cellular senescence. Using an evolution approach, we show that trisomy 8 mitigates EWS-FLI1-induced replication stress through gain of a copy of RAD21. Low-level ectopic expression of RAD21 is sufficient to dampen replication stress and improve proliferation in EWS-FLI1-expressing cells. Conversely, deleting one copy in trisomy 8 cells largely neutralizes the fitness benefit of chromosome 8 gain and reduces tumorgenicity of a Ewing sarcoma cancer cell line in soft agar assays. We propose that RAD21 promotes tumorigenesis through single gene copy gain. Such genes may explain some recurrent aneuploidies in cancer.

Comparative genome-scale analysis of Pichia pastoris variants informs selection of an optimal base strain.

Komagataella phaffii, also known as Pichia pastoris, is a common host for the production of biologics and enzymes, due to fast growth, high productivity, and advancements in host engineering. Several K. phaffii variants are commonly used as interchangeable base strains, which confounds efforts to improve this host. In this study, genomic and transcriptomic analyses of Y-11430 (CBS7435), GS115, X-33, and eight other variants enabled a comparative assessment of the relative fitness of these hosts for recombinant protein expression. Cell wall integrity explained the majority of the variation among strains, impacting transformation efficiency, growth, methanol metabolism, and secretion of heterologous proteins. Y-11430 exhibited the highest activity of genes involved in methanol utilization, up to two-fold higher transcription of heterologous genes, and robust growth. With a more permeable cell wall, X-33 displayed a six-fold higher transformation efficiency and up to 1.2-fold higher titers than Y-11430. X-33 also shared nearly all mutations, and a defective variant of HIS4, with GS115, precluding robust growth. Transferring two beneficial mutations identified in X-33 into Y-11430 resulted in an optimized base strain that provided up to four-fold higher transformation efficiency and three-fold higher protein titers, while retaining robust growth. The approach employed here to assess unique banked variants in a species and then transfer key beneficial variants into a base strain should also facilitate rational assessment of a broad set of other recombinant hosts.

Characterization of the Extracellular Matrix of Normal and Diseased Tissues Using Proteomics.

The extracellular matrix (ECM) is a complex meshwork of insoluble fibrillar proteins and signaling factors interacting together to provide architectural and instructional cues to the surrounding cells. Alterations in ECM organization or composition and excessive ECM deposition have been observed in diseases such as fibrosis, cardiovascular diseases, and cancer. We provide here optimized protocols to solubilize ECM proteins from normal or tumor tissues, digest the proteins into peptides, analyze ECM peptides by mass spectrometry, and interpret the mass spectrometric data. In addition, we present here two novel R-script-based web tools allowing rapid annotation and relative quantification of ECM proteins, peptides, and intensity/abundance in mass spectrometric data output files. We illustrate this protocol with ECMs obtained from two pairs of tissues, which differ in ECM content and cellularity: triple-negative breast cancer and adjacent mammary tissue, and omental metastasis from high-grade serous ovarian cancer and normal omentum. The complete proteomics data set generated in this study has been deposited to the public repository ProteomeXchange with the data set identifier: PXD005554.

Monitoring Error Rates In Illumina Sequencing.

Guaranteeing high-quality next-generation sequencing data in a rapidly changing environment is an ongoing challenge. The introduction of the Illumina NextSeq 500 and the depreciation of specific metrics from Illumina's Sequencing Analysis Viewer (SAV; Illumina, San Diego, CA, USA) have made it more difficult to determine directly the baseline error rate of sequencing runs. To improve our ability to measure base quality, we have created an open-source tool to construct the Percent Perfect Reads (PPR) plot, previously provided by the Illumina sequencers. The PPR program is compatible with HiSeq 2000/2500, MiSeq, and NextSeq 500 instruments and provides an alternative to Illumina's quality value (Q) scores for determining run quality. Whereas Q scores are representative of run quality, they are often overestimated and are sourced from different look-up tables for each platform. The PPR's unique capabilities as a cross-instrument comparison device, as a troubleshooting tool, and as a tool for monitoring instrument performance can provide an increase in clarity over SAV metrics that is often crucial for maintaining instrument health. These capabilities are highlighted.

Comparative genomics and transcriptomics of Pichia pastoris.

BACKGROUND: Pichia pastoris has emerged as an important alternative host for producing recombinant biopharmaceuticals, owing to its high cultivation density, low host cell protein burden, and the development of strains with humanized glycosylation. Despite its demonstrated utility, relatively little strain engineering has been performed to improve Pichia, due in part to the limited number and inconsistent frameworks of reported genomes and transcriptomes. Furthermore, the co-mingling of genomic, transcriptomic and fermentation data collected about Komagataella pastoris and Komagataella phaffii, the two strains co-branded as Pichia, has generated confusion about host performance for these genetically distinct species. Generation of comparative high-quality genomes and transcriptomes will enable meaningful comparisons between the organisms, and potentially inform distinct biotechnological utilies for each species. RESULTS: Here, we present a comprehensive and standardized comparative analysis of the genomic features of the three most commonly used strains comprising the tradename Pichia: K. pastoris wild-type, K. phaffii wild-type, and K. phaffii GS115. We used a combination of long-read (PacBio) and short-read (Illumina) sequencing technologies to achieve over 1000X coverage of each genome. Construction of individual genomes was then performed using as few as seven individual contigs to create gap-free assemblies. We found substantial syntenic rearrangements between the species and characterized a linear plasmid present in K. phaffii. Comparative analyses between K. phaffii genomes enabled the characterization of the mutational landscape of the GS115 strain. We identified and examined 35 non-synonomous coding mutations present in GS115, many of which are likely to impact strain performance. Additionally, we investigated transcriptomic profiles of gene expression for both species during cultivation on various carbon sources. We observed that the most highly transcribed genes in both organisms were consistently highly expressed in all three carbon sources examined. We also observed selective expression of certain genes in each carbon source, including many sequences not previously reported as promoters for expression of heterologous proteins in yeasts. CONCLUSIONS: Our studies establish a foundation for understanding critical relationships between genome structure, cultivation conditions and gene expression. The resources we report here will inform and facilitate rational, organism-wide strain engineering for improved utility as a host for protein production.

RNA polymerase II promoter-proximal pausing in mammalian long non-coding genes.

Mammalian genomes encode a large number of non-coding RNAs (ncRNAs) that greatly exceed mRNA genes. While the physiological and pathological roles of ncRNAs have been increasingly understood, the mechanisms of regulation of ncRNA expression are less clear. Here, our genomic study has shown that a significant number of long non-coding RNAs (lncRNAs, >1000 nucleotides) harbor RNA polymerase II (Pol II) engaged with the transcriptional start site. A pausing and transcriptional elongation factor for protein-coding genes, tripartite motif-containing 28 (TRIM28) regulates the transcription of a subset of lncRNAs in mammalian cells. In addition, the majority of lncRNAs in human and murine cells regulated by Pol II promoter-proximal pausing appear to function in stimulus-inducible biological pathways. Our findings suggest an important role of Pol II pausing for the transcription of mammalian lncRNA genes.

Tumor Cell-Driven Extracellular Matrix Remodeling Drives Haptotaxis during Metastatic Progression.

UNLABELLED: Fibronectin (FN) is a major component of the tumor microenvironment, but its role in promoting metastasis is incompletely understood. Here, we show that FN gradients elicit directional movement of breast cancer cells, in vitro and in vivo Haptotaxis on FN gradients requires direct interaction between α5ß1 integrin and MENA, an actin regulator, and involves increases in focal complex signaling and tumor cell-mediated extracellular matrix (ECM) remodeling. Compared with MENA, higher levels of the prometastatic MENA(INV) isoform associate with α5, which enables 3-D haptotaxis of tumor cells toward the high FN concentrations typically present in perivascular space and in the periphery of breast tumor tissue. MENA(INV) and FN levels were correlated in two breast cancer cohorts, and high levels of MENA(INV) were significantly associated with increased tumor recurrence as well as decreased patient survival. Our results identify a novel tumor cell-intrinsic mechanism that promotes metastasis through ECM remodeling and ECM-guided directional migration. SIGNIFICANCE: Here, we provide new insight into how tumor cell:ECM interactions generate signals and structures that promote directed tumor cell migration, a critical component of metastasis. Our results identify a tumor cell-intrinsic mechanism driven by the actin regulatory protein MENA that promotes ECM remodeling and haptotaxis along FN gradients. Cancer Discov; 6(5); 516-31. ©2016 AACR.See related commentary by Santiago-Medina and Yang, p. 474This article is highlighted in the In This Issue feature, p. 461.

Genetic variations in EGFR and ERBB4 increase susceptibility to cervical cancer.

OBJECTIVES: Inherited genetic variability contributes to susceptibility to cervical cancer. We investigated the association of single nucleotide polymorphisms (SNPs) in the human epidermal growth factor receptor (ERBB) family with cervical cancer. METHODS: We used the transmission disequilibrium test (TDT) to look for excessive transmission of tag single nucleotide polymorphisms (tSNPs) in ERBB family members EGFR, ERBB2, ERBB3, and ERBB4 in a large sample of women with invasive and in situ cervical cancer and their biological parents (628 trios). The study used a discovery set of trios (244) analyzed by Illumina GoldenGate in which SNPs reaching a P<.05 were re-tested by TaqMan in the combined set of 628. We also explored collaborative effects of different ERBB alleles. RESULTS: Based on single SNP TDT tests we identified 16 significant SNPs in the discover stage and six of 14 SNPs that could be assayed by TaqMan were significantly overtransmitted in women with cervical cancer in the combined replication set. Four SNPs were located in intron 1 of EGFR and two SNPs in intron 24 of ERBB4. The EGFR variants are located near multiple enhancers, silencers, and the previously identified functional common polymorphisms in intron 1. CONCLUSIONS: Our data provide evidence for the involvement of intron 1 EGFR variants and intron 24 ERBB4 variants in modulating risk for the development of in situ and invasive cervical cancer. These variants should be examined in additional populations and functional studies would be needed to confirm this hypothesis.

Exome sequencing reveals DNAJB6 mutations in dominantly-inherited myopathy.

OBJECTIVE: To identify the causative gene in an autosomal dominant limb-girdle muscular dystrophy (LGMD) with skeletal muscle vacuoles. METHODS: Exome sequencing was used to identify candidate mutations in the studied pedigree. Genome-wide linkage was then used to narrow the list of candidates to a single disease-associated mutation. Additional pedigrees with dominant or sporadic myopathy were screened for mutations in the same gene (DNAJB6) using exome sequencing. Skeletal muscle from affected patients was evaluated with histochemistry and immunohistochemical stains for dystrophy-related proteins, SMI-31, TDP43, and DNAJB6. RESULTS: Exome analysis in 3 affected individuals from a family with dominant LGMD and vacuolar pathology identified novel candidate mutations in 22 genes. Linkage analysis excluded all variants except a Phe93Leu mutation in the G/F domain of the DNAJB6 gene, which resides within the LGMD locus at 7q36. Analysis of exome sequencing data from other pedigrees with dominant myopathy identified a second G/F domain mutation (Pro96Arg) in DNAJB6. Affected muscle showed mild dystrophic changes, vacuoles, and abnormal aggregation of proteins, including TDP-43 and DNAJB6 itself. INTERPRETATION: Mutations within the G/F domain of DNAJB6 are a novel cause of dominantly-inherited myopathy. DNAJB6 is a member of the HSP40/DNAJ family of molecular co-chaperones tasked with protecting client proteins from irreversible aggregation during protein synthesis or during times of cellular stress. The abnormal accumulation of several proteins in patient muscle, including DNAJB6 itself, suggest that DNAJB6 function is compromised by the identified G/F domain mutations.

HLA-Cw group 1 ligands for KIR increase susceptibility to invasive cervical cancer.

Inherited genetic polymorphisms within immune response genes have been shown to associate with risk of invasive cervical cancer (ICC) and its immediate precursor, cervical intraepithelial neoplasia grade 3. Here, we used the transmission/disequilibrium test to detect disease-liability alleles and investigate haplotype transmission of KIR and HLA class I polymorphisms in a large family-based population of women with cervical cancer and their biological parents (359 trios). The effect of distinct human papillomavirus types was also explored. HLA-Cw group 1 (HLA-Cw alleles with asparagine at position 80), which serves as ligand for certain killer immunoglobulin-like receptors (KIR), was significantly overtransmitted in women with ICC (P = 0.04), and particularly in the subgroup of women infected with high risk HPV16 or 18 subtypes (P = 0.008). These data support the involvement of the HLA-C locus in modulating the risk of cervical neoplasia perhaps through its function as ligands for KIR, but functional studies are essential to confirm this hypothesis.

Loss of Id2 potentiates the tumorigenic effect of Rb inactivation in a mouse model of retinoblastoma.

PURPOSE: In some cancers, the oncogenic consequences of inactivating the retinoblastoma protein (Rb) appear to be mediated by unrestrained activity of the inhibitor of DNA binding protein Id2. The role of Id2 has not yet been investigated in the prototype cancer Rb-defective cancer, retinoblastoma itself. This study investigated whether loss of Id2 modified the effects of Rb inactivation in a mouse model of retinoblastoma. METHODS: Id2 was analyzed in cultured cells using qPCR, Western blot, and colony formation assays. LH beta-Tag transgenic mice were crossed with Id2 heterozygotes to obtain mice with all three Id2 genotypes. Intraocular tumors were assessed for size, degree of differentiation, mitotic index, and tumor vascular density at 15 weeks of age. RESULTS: Retinoblastoma cell lines expressed low levels of Id2 mRNA and protein. Depletion of Id2 in Rb-inactivated cells increased clonogenic activity. Id2-deficient tumors in vivo were significantly larger, less differentiated, and more vascularized than Id2-wild-type tumors (P = 0.02, P = 0.01, P = 0.0001, respectively). There was a dosage effect for loss of each Id2 allele with respect to differentiation and vascular density. CONCLUSIONS: Id2 suppresses rather than promotes tumor progression in this mouse model of retinoblastoma. Id2 can act as either an oncogene or a tumor suppressor depending on context.

TP53, MDM2, NQO1, and susceptibility to cervical cancer.

Host genetic variability modifies the risk of cervical cancer in women infected with oncogenic human papillomavirus (HPV). Studies have reported an association of the TP53 codon 72 arginine and cervical cancer, but the results are inconsistent. We examined the association of this single nucleotide polymorphism (SNP) in women with cervical cancer and cervical intraepithelial neoplasia grade 3, using a family-based association test. We further explored SNPs in two genes that regulate p53 stability: MDM2 (SNP309) and NQO1 (SNP609, SNP465). We also examined the relationship between host genotype and tumor HPV type. We genotyped 577 patients and their biological parents and/or siblings, using PCR-RFLP or Taqman assays. HPVs were typed by sequence-based methods. The transmission/disequilibrium test was used to detect disease-susceptibility alleles. The arginine peptide of TP53 codon 72 was overtransmitted in Caucasian families (P = 0.043), and the significance of this finding was enhanced in a subgroup of women infected with HPV16- and/or HPV18-related HPVs (P = 0.026). Allele C of NQO1 SNP609 was also overtransmitted in all cases (P = 0.026). We found no association between MDM2 SNP309 or NQO1 SNP465 and cervical cancer. Our results indicate that functional polymorphisms in TP53 codon 72 and NQO1 SNP609 associate with the risk of cervical cancer especially in women infected with type 16- and/or type 18-related HPVs.

Longitudinal trends in the association of metabolic syndrome with 550 k single-nucleotide polymorphisms in the Framingham Heart Study.

We investigated the association of metabolic syndrome (MetS) with a 500 k and a 50 k single-nucleotide polymorphism (SNP) gene chip in the Framingham Heart Study. We cross-sectionally evaluated the MetS longitudinal trends. Data analyzed were from the Offspring Cohort (four exams: first (n = 2,441), third (n = 2,185), fifth (n = 2,308), and seventh (n = 2,328)) and the Generation 3 Cohort (one exam: the first exam (n = 3,997)). The prevalence of MetS was determined using the National Cholesterol Education Program Adult Treatment Panel III diagnostic criteria, modified with a newly developed correction for medication use. The association test between an SNP and MetS was performed with a generalized estimating equations method under the additive genetic model. Multiple-testing corrections were also performed. The prevalence of MetS in the offspring cohort increased from one visit to the next, and reached the highest point by the seventh exam comparable with the prevalence reported for the general US population. The pattern of the MetS prevalence over time also reflected itself in the association tests, in which the highest significances were seen in the fifth and seventh exams. The association tests showed that SNPs within genes PRDM16, CETP, PTHB1, PAPPA, and FBN3, and also some SNPs not in genes were significant or close to significance at the genome-wide thresholds. These findings are important in terms of eventually identifying with the causal loci for MetS.

A framework for analyzing both linkage and association: an analysis of Genetic Analysis Workshop 16 simulated data.

We examine a Bayesian Markov-chain Monte Carlo framework for simultaneous segregation and linkage analysis in the simulated single-nucleotide polymorphism data provided for Genetic Analysis Workshop 16. We conducted linkage only, linkage and association, and association only tests under this framework. We also compared these results with variance-component linkage analysis and regression analyses. The results indicate that the method shows some promise, but finding genes that have very small (<0.1%) contributions to trait variance may require additional sources of information. All methods examined fared poorly for the smallest in the simulated "polygene" range (h2 of 0.0015 to 0.0002).